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Abstract 

Drought is one of the most severe stresses, endangering crop yields worldwide. In order to select drought tolerant 
genotypes, access to exotic germplasm and efficient phenotyping protocols are needed. In this study the high-throughput 
phenotyping platform "The Plant Accelerator", Adelaide, Australia, was used to screen a set of 47 juvenile (six week old) wild 
barley introgression lines (S42ILs) for drought stress responses. The kinetics of growth development was evaluated under 
early drought stress and well watered treatments. High correlation {r = 0.98) between image based biomass estimates and 
actual biomass was demonstrated, and the suitability of the system to accurately and non-destructively estimate biomass 
was validated. Subsequently, quantitative trait loci (QTL) were located, which contributed to the genetic control of growth 
under drought stress. In total, 44 QTL for eleven out of 14 investigated traits were mapped, which for example controlled 
growth rate and water use efficiency. The correspondence of those QTL with QTL previously identified in field trials is 
shown. For instance, six out of eight QTL controlling plant height were also found in previous field and glasshouse studies 
with the same introgression lines. This indicates that phenotyping juvenile plants may assist in predicting adult plant 
performance. In addition, favorable wild barley alleles for growth and biomass parameters were detected, for instance, a 
QTL that increased biomass by approximately 36%. In particular, introgression line S42IL-121 revealed improved growth 
under drought stress compared to the control Scarlett. The introgression line showed a similar behavior in previous field 
experiments, indicating that S42IL-121 may be an attractive donor for breeding of drought tolerant barley cultivars. 
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Introduction 

Barley [Hordeum vulgare ssp. vulgare, hereafter abbreviated with 
Hv) is ranked fourth among the worldwide production of cereals. 
Due to its multipurpose use as animal feed, human food and 
substrate for malting it is one of the most important cereals world- 
wide [1]. Barley is known to be relatively tolerant to abiotic stresses 
among the major cereal crops and, thus, is often grown in more 
marginal sites [2]. However, the process of genetic erosion has 
been under way in barley since its domestication some 10,000 
years ago and, in particular, since the advent of intensive modern 
elite breeding during the last century [3]. As a result of this 
process, diverse landraces have been replaced by modern elite 
cultivars with a much narrower gene pool. Therefore there is 
limited genetic diversity remaining in the elite barley gene pool for 
abiotic and biotic stress tolerance. The current loss of genetic 
variation in the elite gene pool tends to limit the breeding success 
of improved cultivars [4]. To overcome this problem several 
authors, e.g. Zamir [5], proposed to use wild relatives of crop 
species as donors of exotic germplasm to enhance elite varieties. 
Tanksley and Nelson [6] proposed the method of "advanced 



backcross quantitative trait loci analysis" (AB-QTL) to introduce 
exotic genes into modern crop varieties. The method combines 
QTL detection and the introduction of favorable exotic alleles 
from a wild donor parent. Lines produced by advanced 
backcrosses ideally contain only one single introgxession from 
the wild parent and are then referred to as introgression lines (ILs). 
This is achieved by several rounds of backcrossing to the recurrent 
parent and marker assisted selection (MAS). A set of ILs ideally 
represents the whole genome of a wild donor plant in the genetic 
background of a single elite variety [5] . 

Pillen et al. [7] published the first AB-QTL study in barley. Von 
KorfF et al. [8] developed a BC2DH population from a cross 
between the German spring barley cultivar Scarlett [Hv) and the 
Israeli wild barley accession ISR42-8 [Hordeum vulgare ssp. 
spontaneum, hereafter abbreviated with Hsp). The lines of this S42 
population were used in several AB-QTL studies to identify QTL 
for yield, pathogen resistance and malting quality traits [9-13]. 

Schmalenbach et al. [14] used the S42 population to develop 59 
ILs (S42ILs) by a further round of backcrossing with the recurrent 
parent Scarlett and subsequent selling and MAS. Each of the 
S42ILs contains a single or a small number of Hsp introgressions. 
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Several QTL studies were conducted to verify QTL from AB- 
QTL studies and to identify new QTL for pathogen resistance, 
yield and quality parameters [14—16]. Naz et al. [17] studied root 
architecture of S42ILs and detected QTL for root dry weight and 
root volume. Later on, the S42IL population was extended to 73 
lines and the lines were genotyped with a 1,536-SNP lUumina 
BOPAl set [18]. Six hundred thirty-six informative SNPs and 
their known map order [19] allowed the precise localization of the 
Hsp introgressions. The S42IL set represents 87.3% of the wild 
barley donor genome. Moreover, Schmalenbach et al. [18] 
developed segregating high-resolution mapping populations 
(S42IL-HRs) for 70 S42ILs. Those lines are readily available to 
facilitate fine mapping and, ultimately, cloning of QTL. 

Drought is one of the main factors limiting yield worldwide 
[20] . Due to climate change extreme weather events are predicted 
to occur more frequently and an altered pattern of drought 
occurrence is expected [21]. Therefore maintaining plant growth 
and yield under drought remains a major objective for plant 
breeding [22]. Many studies have been conducted on the impacts 
of terminal drought stress in cereals, while impacts of drought 
stress at early developmental stages are less well investigated [23]. 
However, several authors comment that yield may be enhanced by 
improved early vigor and a rapid development of maximum leaf 
area [24,25]. Lopez-Castaiieda and Richards [26] reported that 
on average barley has a higher yield in water limited environments 
compared to wheat, triticale, and oat. As part of a possible 
explanation, they pinpointed the faster and more vigorous growth 
of barley during vegetative development. Variation in this trait is, 
therefore, likely to be in direct relation to drought stress tolerance 
and yield. 

Conventional methods to determine biomass and measure 
growth are time-consuming and labor intensive. Often they 
involve destructive harvest of plants and therefore make repeated 
measurements on the same plant impossible. New developments in 
plant imaging technologies allow the estimation of biomass and 
growth parameters as a non-destructive and rapid alternative to 
more traditional methods [27,28]. New phenotyping facilities 
enable automated imaging of plants. Several types of plant images 
can be taken, e.g. with infrared, near infrared, fluorescent and 
visible light. Scanning with infrared light gives information on 
plant or leaf temperature, while near infrared imaging sheds light 
on the plant water content and fluorescent pictures enable 
conclusions on plant health status. High resolution color pictures 
(RGB pictures), taken from the top and two side views are used to 
determine the projected shoot area of the plant. The projected 
shoot area serves as a measure for biomass. Hence, from RGB 
images taken at several time points, growth curves as well as 
growth rates can be calculated. 

In fully automated greenhouses plants can be delivered via 
conveyor belts to watering, weighing and imaging stations. In 
these high-throughput phenotyping facilities several hundred 
individual plants can be imaged per day in a fully automated 
manner. High-throughput phenotyping facilities of this type are 
currently in use in various research institutes (e.g. The Plant 
Accelerator, Adelaide, Australia; CropDesign, Gent, Belgium; IPK 
Gatersleben, Germany, PhenoArch, MontpeUier, France). 

Such phenotyping facilities are ideal to combine controlled 
irrigation and phenotyping protocols [29] . A first application was 
given by Rajendran et al. [30] who used a manual imaging system 
(LemnaTecScanalyzer3D, Wuerselen, Germany) to screen Triticum 
momcoccum accessions for salinity tolerance. They developed high- 
throughput quantification assays to distinguish sodium exclusion, 
sodium tissue tolerance and osmotic tolerance as the strategies 
plants use to establish salinity tolerance. 



In this report, "The Plant Accelerator" was used to screen 
growth of wild barley ILs under well watered and drought 
treatments during vegetative growth. The aims were (1) to identify 
wild barley derived QTL within the set of S42ILs that control 
drought stress responses and (2) to test the use of non-destructive 
high-throughput imaging to measure vegetative stage drought 
response in barley. 

We could show that high-throughput imaging provides accurate 
estimates for biomass development over time. Moreover several 
drought related QTLs were identified and genotypes detected that 
may be beneficial in future breeding programs. 

Materials and Methods 

Plant Material 

Forty-seven wild barley ILs of the S42IL library and the 
recipient parent Scarlett were selected for the experiment. The 
S42ILs are derived from a cross between the German malting 
barley variety Scarlett and the Israeli wild barley accession ISR42- 
8. The 47 ILs possess few H.sp chromosome segments and were 
selected based on SSR and SNP genotyping to represent a large 
portion, 87.3%, of the ISR42-8 genome [18]. Repeated 
backcrossing and MAS are described in Schmalenbach et al. [14]. 

Glass House Cultivation 

Two drought stress experiments, with duration of six weeks 
each, were conducted between end of March and mid of July 201 1 
in The Plant Accelerator greenhouse facilities in Adelaide, 
Australia (34°58' 16.18"S; 138°38'23.88"E). Forty-eight barley 
genotypes were grown under a well watered and stress treatments 
with three replicates per genotype and treatment. Each experi- 
ment was designed in three randomized blocks. Control and stress 
treatments of each genotype were placed next to each other (Fig. 1, 
Table SI). 

Single plants were grown in 2.5 L plastic containers with 2.1 kg 
of soil (50% UC Davis soil mix, 35% Coco-peat, 15% clay-loam). 
Three seeds per pot were directly sown into the soil and after 
germination thinned out, leaving one plant per pot. Plants were 
pre-grown for two weeks in a regular greenhouse and watering was 
performed manually to allow optimal germination and seedling 
estabhshment. Subsequentiy, the pots were transferred to the 
"smart house" where each pot was placed onto a cart on a 
conveyor belt and the two treatments were applied. Every second 
day, pots were weighed and watered automatically to 22% 
gravimetric water content for the well watered treatment and 1 5 % 
for the stress treatment (Fig. 2). Based on the experience of the first 




Figure 1. View of experiment 1 with five-week old barley S42IL 
plants growing in The Plant Accelerator. 

doi:10.1371/journal.pone.0097047.g001 
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Figure 2. Barley plants at the weighing and watering unit after 
leaving the imaging station. 

doi:1 0.1 371 /journal.pone.0097047.g002 

experiment we adjusted the drought stress in the second 
experiment to 12% gravimetric water content to shghtly increase 
drought effects. The experiments were carried out under natural 
lighting with the temperature in the greenhouse kept at a range 
between 15°C (night) and 22°C (day). 

Phenotyping 

With the onset of the stress treatment imaging of the plants 
started. Plant images were captured using a LemnaTec 3D 
Scanalyzer (LemnaTec, GmbH, Wuerselen, Germany). Every day, 
three RGB pictures (2056x2454 pixels) were taken of each barley 
plant, one top view image and two side view images with a 90° 
horizontal rotation. After background-foreground separation was 
applied to separate the plant tissue area from the background, 
pixel numbers per plant were counted and the pixel sum of the 
three pictures per plant was taken to define the projected shoot 
area. The shoot area measured over time was used to draw growth 
curves. For each growth curve, curve fitting with a 6* order 
polynomial was conducted to adjust for possible missing data 
points and absolute growth rate [dA/ dt] and relative growth rate 
[(dA/dt)/A] were calculated. For each of the three curves the 
integral was determined and used as a trait in the statistical 
analysis. Moreover, six further traits were extracted from the 
images; caliper length, height, color (as hue angle in the HSI color 
scheme) and the two parameters shoot area top view and convex 
huU area to calculate compactness of each plant. At the end of the 
experiment, barley plants were harvested and above ground 
biomass, tiller number (TIL), and plant height (HEI) were 



determined. Fresh biomass was weighed and, subsequendy, oven 
dried to constant weight to determine dry biomass. Water use 
efficiency (WUE) was calculated by dividing dry biomass at the 
end of the experiment by the total amount of water added during 
the four weeks in the "smart house" [mg/g water]. Specific plant 
weight (SPW) was calculated from the dry weight and the 
maximum projected shoot area at the end of the experiment. In 
addition, simple stress indices (SSI) were calculated as follows: 
SSI = Ts/Tc, where Ts and Tc are the average trait performances 
of an IL under stress and control conditions, respectively. An 
overview of trait definitions is given in Table 1 . 

Genotyping 

The S42ILs were genotyped with the 1 ,536-SNP bariey BOPAl 
set [19] of the lUumina GoldenGate assay [18]. Six hundred and 
thirty-six out of the tested 1,536 SNPs were polymorphic and used 
to characterize the extent of exotic Hsp introgressions in each 
S42IL (see Fig. SI). 

Statistical Analysis 

Statistical analyses were performed with SAS Enterprise 
Guide 4.2. [31]. Descriptive statistical parameters 
(Table S2) were calculated with procedure MEANS. 
HeritabHities across treatments were calculated as 
h^=VGl[VG+VGTlt-\- VGE/e+ VGET/et+ VR/etr] , and withm 
treatments as: = Kg/[Kg + KG£/e+ K^/er]. The terms Vq, 
Vgt, Vge, Vget and Vr represent the genotypic, genotype x 
treatment, genotype Xenvironment, genotype x environment x 
treatment, and error variance components, respectively, calculated 
with procedure VARCOMP [31]. The terms t, e, and r indicate 
the number of treatments, experiments and replicates, respective- 
ly. Pearson correlation coefficients between traits were calculated 
with means across treatments, blocks and experiments and within 
drought stressed and control treatments, respectively, using the 
procedure CORR. 

Analysis of variance was carried out with the procedure 
MIXED using model I to test for genotype main effects across 
treatments and experiments. 

Model I: 

Yijki = + L, + Tj + Ek + LxT,i + Lx + B{Ex Tkj) , + Eyu 

and model II for genotype effects across experiments but within 
a single treatment. 
Model II: 

Yu = n + L, + Ek + Lx E,k + B{Ek)i + Em 

Where jX is the general mean, i, is the fixed effect of the ith line, Tj 
is the frxed effect of the jth treatment, E^ is the fixed effect of the 
bh experiment, L x Tij is the fixed interaction between the Jth Une 
and thejth treatment, LxEn^ is the fixed interaction between the 
rth line and the Ath experiment, BfExT/^Ji is the random effect of 
the /th block nested in the interaction between kth experiment and 
jth treatment, BfE^ji is the random effect of the /th block nested in 
the ^h experiment and and £,h are the error of I'j^a a^nd T^^i, 
respectively. 

Following the mixed model analysis a Dunnett test was 
conducted where least square means (LSMEANS) of each IL 
were compared to the control Scarlett. In case an IL revealed a 
significant (P<0.05) deviation in trait performance from Scarlett, 
as main effect and/or as line x treatment interaction, a line x trait 
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Table 1. List of evaluated traits. 





Trait 


Abbreviation 


Unit 


Method of measurement 


Imaging parameters 


Shoot area integral^ 


SAI 


kPix^ 


Calculated from pixel sum of three images per plant per day; A 


Absolute growth rate integral 


AGRI 


kPix/d 


Calculated from pixel sum of three images per plant per day; dA/dt 


Relative growth rate integral 


RGRI 


d-^ 


Calculated from pixel sum of three images per plant per day; (dA/dt)/A 


Height integral 


HEN 


kPix 


Max. distance from bottom to top of plant 


Caliper length integral 


CALI 


kPix 


Max. distance between two points on the object boundary, top view image 


Hull area integral 


HULI 


kPix 


Smallest geometrical object without concave parts that covers whole plant, top view image 


Shoot area top view integral 


SATVI 


kPix 


Pixel number 


Plant hue integral 


HUEI 




Average hue value calculated from all pixels per plant and day 


Harvest parameters 


Tiller number 


TIL 




Number of tillers per pot 


Height 


HEI 


cm 


Plant height measured from bottom to leaf tip 


Biomass dry 


BMD 


g 


Weight of oven dried biomass per pot 


Ir) dices 


Water use efficiency 


WUE 


mg/g 
water 


Harvested biomass per plant/total amount of irrigation water 


Specific plant weight 


SPW 


mg/kPix 


Harvested biomass per plant/pixel number per plant at end of experiment 


Compactness integral 


COMI 




SATV/HUL per plant per day 


Simple stress index 


SSI 




Trait performance under stress/trait performance under control treatment 



^Integral: calculated for length of entire experiment, respectively, 
"kilo Pixel. 

doi:l 0.1 371/journal.pone.0097047.t001 



association was assumed and tlie presence of a QTL was accepted. 
If several ILs witti overlapping introgressions showed a similar 
effect, it was assumed that the ILs contained the same QTL. We 
consider this QTL as the most likely location of the effect and, 
thus, defme a minimum number of QTL needed to explain all 
identified trait effects. The relative performance (RP) of an IL was 
calculated as RP (IL) = [LSMEANS (IL) - LSMEANS (Scariett)] 
xlOO/LSMEANS (Scariett), where LSMEANS were calculated 
with model I across treatments, experiments and blocks or with 
model II across experiments and blocks, separately for each 
treatment. The detection of significant line by trait associations 
was conducted for every trait revealing a heritabUity with h^>0% 
across treatments or within the two watering treatments, 
respectively. 

Results 

Trait Performance of S42ILs 

For most traits, means were higher under well watered 
treatment than under drought stress (Fig. 3 and Fig. 4, Table 
S2), e.g. 2.2 g of biomass dry (BMD) vs. 0.9 g. There were, 
however, four exceptions. Compactness integral (COMI) was 
higher under drought treatment than under well watered 
treatment. The same was true for plant hue integral (HUEI), 
SPW, and WTJE but differences were marginal. 

Coefficients of variation (CV) differed strongly between traits 
(1.5 to 72.1%). Highest CV was calculated for BMD across 
treatments (72.1%). The lowest CVs were determined for HUEI, 
varying from 1.5 to 1.6% for the different treatments and the SSI. 
CV was generally higher under drought than under well watered 
treatment. The four exceptions were WUE, SPW and height 
integral (HEII) where CV under well watered treatment was 
higher and HUEI where CV was the same under both treatments. 



Heritability was generally higher under well watered treatment. 
HEI, HEII and WUE were exceptions and showed higher 
heritabUities under drought treatment. Highest heritabilities were 
found for HEI, HEII and caliper integral (CALI) (between 46.2 
and 76.8%). Low heritabilities were determined for WUE and 
HUEI under drought treatment and across treatments, as well as 
for BMD under drought treatment (15%). Most of the SSI showed 
heritabilities equal to 0. SSI (HEII), SSI (HUEI), SSI (SPW), SSI 
(WUE) revealed heritabilities between 4.7 and 14.5%. 

Trait Correlations 

Highest correlations were found among measured traits and 
among stress indices (Table S3). However, correlations between 
measured traits and stress indices were low. The measured traits 
showed the highest correlations between shoot area integral (SAI) 
and absolute growth rate integral (AGRI) (r = 0.99), BMD and SAI 
(Fig. 4), BMD and AGRI, shoot area top view integral (SATVI) 
and SAI and SATVI and AGRI (r = 0.98). Most correlations were 
positive and statistically significant. HEI showed negative corre- 
lation with relative growth rate integral (RGRI) and TIL, but 
values were not statistically significant. COMI showed negative 
correlations with all traits but HEII. Among stress indices highest 
correlations were found for SSI (SAI), SSI (AGRI) and SSI (SAI), 
as weU as between SSI (WUE) and SSI (SPW) (with r>0.93). 
Interestingly, WUE, SPW and BMD showed only low correlations 
(r<0.43). Looking at the simple stress index, however, correlations 
between SSI (WUE), SSI (SPW) and SSI (BMD) were very high 
(r = 0.85 to 0.96). Autocorrelations between drought and well 
watered treatments of a single trait were high (r>0.61) for most 
traits. HEI with r = 0.85 had the highest correlation between 
treatments. RGRI showed a low but still significant correlation 
between the treatments with r = 0.33. Autocorrelations for SPW 
and WUE were not significant. 
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Figure 3. Development of shoot area of S42IL-121, Scarlett and S42IL-1 17 under well watered (solid line) and drought (dashed line) 
treatment, respectively. 

doi:1 0.1 371 /journal.pone.0097047.g003 



Mixed Model Analysis of Variance 

The mixed model analysis, including fixed line, treatment, and 
experiment effects (i.e. model I) revealed significant (P<0.05) line 
effects for aU investigated traits (Table S4). Treatment had a clear 
impact on trait performance. For all traits, except WUE and SPW 
the effect was significant. Line by treatment interaction effects 
were not significant for any of the measured traits. The 
experiments had a significant effect on trait performance of all 
traits except leaf color measured as plant hue integral (HUEI). And 
line by experiment interaction was significant for all traits but 
RGRI, COMI, and HEII. In the mixed model analyses for single 
treatments including fixed line, and experiment effects (i.e. model 
II), line had a significant effect on trait performance for aU traits 
but HUEI, RGRI, and WUE under well watered and HUEI and 
SPW under drought conditions (Tables S,5 and S6). Also simple 



stress indices were analyzed with model II (Table S7). The line 
effect was not significant for any of the simple stress indices. 

QTL Detection 

QTL were only determined for traits with heritability greater 
than 0. The Dunnett tests revealed, in total, 63 line effects for 
eleven out of 14 traits. These effects were detected either across 
treatments (39), within the drought treatment (15) or within the 
well watered treatment (9). Several of the measured effects were 
consistent between the different treatments. Thus, these line effects 
were summarized to a minimum of 44 QTL (Table 2 and Fig. 5). 
No QTL were identified for RGRI, SPW, and HUEI. Between 
two and nine QTL were identified for the traits AGRI, BMD, 
CALI, COMI, HEI, HEII, hull area integral (HULI), SAI, 



(r = 0.98) 



1* 



watered 



J^ro 



drought 



1000 2000 3000 

Shoot Area in kiloPixel 



4000 



Figure 4. Correlation between biomass and shoot area integral under drought (blue dots) and well watered (red dots) treatment, 
respectively. 

doi:1 0.1 371/journal.pone.0097047.g004 
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Figure 5. QTL map with indication of S42IL introgressions (Schmalenbach et al. 2011). SNP positions (in cM) are based on Close et al. 
(2009). QTL are placed right to the S42IL, indicated by trait abbreviations (see Table 1). The sign indicates an increasing (+) or a decreasing (— ) Hsp 
effect. 

doi:1 0.1 371 /journal.pone.0097047.g005 
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SATVI, TIL, and WUE. In the foUowing, the QTL are presented 
for each trait separately. 

Absolute Growth Rate Integral (AGRI) 

Three QTL were identified for AGRI. The QTL are located on 
chromosomes 3H, 4H and 6H and the Hsp allele in aU three QTL 
reduced the trait performance. Across treatments the Hsp alleles 
reduced the integral of the absolute growth rate by approximately 
30%. Under drought conditions the Hsp allele of QTL 
QAgri.S42IL-4H reduced the trait performance by almost 40%. 

Biomass Dry (BMD) 

For biomass four QTL were identified across treatments. The 
Hsp alleles at QTL QBmd.S42IL-3H and QBmd.S42IL-6H on 
chromosomes 3H and 6H reduced biomass by approximately 
33%. The two QTL QBmd.S42IL-4H and QBmd.S42IL-4Hb on 
chromosome 4H showed contrary elfects. While at the first QTL 
the Hsp allele reduced biomass by 40.3%, at the second QTL it 
increased it by 36.0%. 

Caliper Length Integral (CALI) 

Six QTL were detected for caliper length on chromosome IH, 
2H, 3H, 4H and 6H. AU QTL were detected across treatments, 
three and two of them also showed effects for drought and well 
watered treatments, respectively. In five c:ases the Hsp allele had 
decreasing effects between 14.8 and 22.fS% across treatments. In 
one case the Hsp allele at QTL QCali.S42IL-4Hb increased 
caliper length by 15.2% compared to Scarlett. 

Compactness Integral (COMI) 

For COMI two QTL were detected on chromosomes 4H and 
6H. The effect of the Hsp allele at QComi.S42IL-4H was observed 
across treatments and under well watered treatment. It increased 
compactness across treatments by 29.8%. In case of QComi.S42IL 
6H the effect was solely observed under well watered treatment 
and lead to an increase of 27.1% by the Hsp allele. 

Height (HEI) 

The highest number of QTL was detected for plant height. 
Eight QTL were detected across treatments. Of those, four also 

showed effects under either one or both of well watered and 
drought treatment. The QTL are located on all chromosomes 
except 5H. In five cases the Hsp alleles reduced plant height (9.2 to 
12.8%). In three cases the Hsp allele increased plant height by 1 1.6 
to 18.7%. 

Height Integral (HEN) 

Five QTL across treatments were detected for HEII on IH, 4H 
and 7H. Two of diose, QHeii.S42IL-lH and QHeii.S42IL-lHb, 
coincided with QTL for manual measurement of HEI. At all 
detected QTL the Hsp alleles reduced HEII by between 13.9 to 
15.6%. 

Hull Area Integral (HULI) 

For HULI six QTL were detected on chromosomes IH, 2H, 

3H, 4H, and 6H. In aU cases the Hsp alleles reduced the hull area 
by between 30.9 to ?>1 A% across treatments. In addition, 
QHuli.S42IL-lH and QHuli.S42IL-2H showed effects under weU 
watered and drought treatments, respectively. 

Shoot Area Integral (SAl) 

Three QTL were found for the integral of the projected shoot 
area. The QTL are located on chromosomes 2H, 4H, and 6H. AU 



three QTL were detected under drought treatment, while, in 
addition, QSai.S42IL-4H and QSai.S42II^6H showed effects 
across treatments. The presence of the Hsp aUele reduced the 
projected shoot area between 30.4 and 37.1% across treatments. 

Shoot Area Top View Integral (SATVI) 

Three QTL were detected for shoot area top view that 
corresponded to the same QTL detected for SAI on chromosomes 
2H, 4H, and 6H. However, only QSatvi.S42IL-2H was detected 
under drought treatment, whUe the two other QTL were detected 
across treatments. In aU cases the Hsp aUeles reduced SATVI 
between 33.6 and 46.2% compared to Scarlett. 

Tiller Number (TIL) 

For tiller number two QTL were identified on chromosomes 3H 
and 4H. The Hsp aUele at QTU.S42IL-3H reduced tUler number 
by 27.6% across treatments. QTil.S42IL-4H was detected across 
treatments and separately within the two treatments. Across both 
treatments the Hsp aUele increased the tiUer number by 42.9%. 

Water Use Efficiency (WUE) 

Two QTL were detected for XWE. The Hsp allele of QTL 
QWue.S42IL-4H and QWue.S42IL-6H reduced WUE by 36.5 
and 40.6% under drought treatment, respectively. 

Relative Growth Rate Integral (RGRI) and Plant Hue 
Integral (HUEI), Specific Plant Weight (SPW), and Simple 
Stress Index (SSI) 

For RGRI, HUEI, and SPW and the SSIs no QTL were 
detected in this study. 

Discussion 

The aim of the study was to validate the use of non-destructive 
high-throughput phenotyping to measure vegetative drought 
response in barley and to identify QTL derived from wUd barley 
that control physiological traits related to drought stress. To the 
authors knowledge this is the first QTL report on drought stress 
that used a high-throughput phenotyping facility. 

Plant growth and the biomass parameters tiller number, plant 
height, and shoot dry weight of 48 barley genotypes were 
investigated under drought and well watered treatments. Two 
week old barley plants were transferred into a high-throughput 
j)lu'iiotyping greenhouse, where stress and control treatments were 
applied automaticaUy. During the foUowing four weeks of 
cultivation, plants were imaged daUy in an automated manner. 
Images were processed and used as a measure for plant height, 
caliper length, biomass and, consequentiy, plant growth, and plant 
color. After a total of six weeks green plants were harvested. TiUer 
number, plant height, and shoot dry weight were determined for 
each plant. Moreover, plant compactness, water use efficiency and 
specific plant weight, as weU as stress indices were calculated. 

The mixed model ANOVA revealed a clear effect of the 
treatment on trait expression. Drought stressed plants had a lower 
growth rate and subsequently produced less biomass (see Fig. 4). 
However, there was no significant line xtreatment interaction, 
indicating that the S42ILs reacted simUar under drought stress and 
weU watered conditions for the investigated traits. This finding is 
also supported by high autocorrelations for investigation of traits 
under drought and control treatments, e.g. 0.77 for biomass (see 
Table S3). 
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QTL Detection 

In this study 44 QTL were identified in 15 ILs for eleven traits. 
In eight cases the Hsp alleles increased the performance of the trait, 
while in 36 cases there was a decreasing effect. This is to be 
expected since wild barley is known to carry many unfavorable 
alleles as well [5]. In six ILs only one QTL was determined, 
predominantly for HEI and HEII. Multiple QTL effects were 
found in nine ILs. For S42IL-1 15,-1 17, -121, and -129 QTL for 
BMD were detected in combination with one or more of the traits 
AGRI, HEI, SAI, and TIL which is in agreement with the high 
correlations found between those traits. SAI is a measure for 
biomass, and AGRI is the first derivative of SAI. Therefore, it was 
expected that ILs show effects for all three traits, simultaneously. 
However, this was not always the case. Thus, increasing the 
number of experiments and replications might be useful to 
increase the power of QTL detection. In the following, the traits or 
trait complexes are discussed separately. 

Absolute Growth Rate Integral, Shoot Area Integral and 
Biomass Dry (AGRI, SAI, BMD) 

The three traits were highly correlated with each other (r = 0.98 
and r = 0.99). Since their relation may be functional, it appears 
likely that a single pleiotropic QTL may control the three traits 
AGRI, SAI and BMD, simultaneously. 

For line S42IL-129 a biomass reduction of 33.9% was observed. 
March et al. (in prep.) found a similar decrease in biomass in that 
line measured under terminal drought stress. This suggests that 
biomass production may be partly controlled by similar genes 
during early and late drought stress occurrence. 

Three QTL were detected for SAI. Two of those, namely 
QSai.S42IL-4H and QSai.S42IL-6H, in line S42II^117 und 
S42IL-129, respectively, were in accordance with BMD QTL. 
The Hsp allele, in both lines caused a decrease in the projected 
shoot area. Due to the high correlation of the traits it can be 
assumed that QTL for biomass correspond to QTL for shoot area. 
The Hsp allele of the QTL QSai.S42IL-2H on chromosome 2H 
caused a decrease of 37% in projected shoot area, compared to 
Scarlett. Von Korff" et al. [10] described a QTL related to biomass 
reduction (Qrv'Ias.S42-2H.a) in the same region in an AB-QTL 
field trial. Since this AB-QTL population is the parent population 
of the S42ILs used in this study, it is likely that the same QTL was 
detected in both the greenhouse and field trials. 

AGRI is directly related to SAI. This might be seen as a reason 
for the detection of similar QTL for AGRI and SAI and, 
consecjuently, BMD. Many QTL studies on growth focus on 
relative growth rate instead of absolute growth rate. In the present 
study RGRI showed only a weak correlation to AGRI (r = 0.35) 
and other biomass parameters. Poorter et al. [33] pointed out that 
in their study QTL for RGR rather co-located with QTL for seed 
mass than with QTL for biomass. This fits well to the weak 
correlations found between biomass parameters and RGRI. 
However, two of the QTL detected for AGRI coincided with 
locations where previous studies mapped QTL for RGR. Yin et al. 
[34] reported a minor effect for relative growth rate associated 
with the denso locus on chromosome 3H in a spring barley 
recombinant inbred lin(; (RIL) population of the cross Prisma x A- 
pex, which may co-localize with QAgri.S42IL-3H. Poorter et al. 
[33] conducted QTL studies in a F2 population derived from a 
cross between two Hsp accessions. They mapped QTL for relative 
growth rate on chromosomes IH, 2H, 5H and a minor QTL on 
6H. The latter one might be in accordance with QAgri.S42IL-6H. 



Tiller Number (TIL) 

Two QTL were detected for the trait tiller number on 
chromosomes 3H and 4H. Wang et al. [35] identified the VKN- 
H'2 gene on chromosome 4H in introgrcssion line S42IL-124. 
Whereas S42IL-124 carries a dominant winter-type allele, Scarlett 
carries the recessive and deleted spring type allele at Vm-H2. 
S42IL-124 showed an increased tiller number compared to 
Scarlett. Since Vrn-H2 is known to have a pleiotropic effect on 
tiller number [36], we assume that this gene explains the 
underlying effect of the QTL. Studies on other populations 
revealed QTL for tiller number on chromosomes 4H as well. In a 
cross between two wild barley accessions Elberse et al. [37] found 
a QTL for tiller number on that chromosome. Baum et al. [38] 
identified a QTL on chromosome 4H where the Hsp allele 
increased the number of tillers and a QTL on chromosome 3H 
where the Hsp allele had a decreasing impact in an ArtaxHsp 41-1 
RIL population. Those effects might correspond to the QTL 
detected in this study. Both QTL occurred irrespective of the 
treatment. Especially QTil.S42IL-4H appears to be a very stable 
QTL. It was detected across and within treatments and was 
detected in several studies under varying conditions, in field studies 
as well as under greenhouse conditions. Moreover, von Korff et al. 
[10] detected QTL for number of ears, which is directiy related to 
tiller number, in the same genomic regions. On 4H the Hsp allele 
increased the number of ears, while on 3H it has a decreasing 
effect. This supports the observation of a stable QTL. 

Height (HEI) and Height Integral (HEII) 

Plant height was determined in two ways. First, height (HEII) 
was modeled from the images taken during four weeks and the 
integral of the height growth curve was calculated. Second, height 
(HEI) was measured manually when plants were harvested after 
six weeks at the end of the experiment. The correlation between 
HEI and HEII was relatively low with r = 0.72, compared to the 
correlation between SAI and BMD with r = 0.98. While SAI 
shows a constant increase over time, HEII shows an overall 
increase, but fluctuation between days may be strong. When a new 
leaf is unfolded the plant grows higher, however, when the leaf 
becomes too heavy and bends down, the height of the plant 
appears to be shorter. At the end of the experiment the length of 
the stretched plant was measured, which is longer than the upright 
standing plant. Nevertheless two coinciding QTL were found 
between HEI and HEII on chromosome IH. 

Six out of eight QTL were already identified in previous field 
studies with the S42 population and/or the S42ILs, exhibiting 
similar effects of the same direction in all three studies. QTL 
QHei.S42IL-3Hb was already detected in von Korff et al. [10] 
and March et al. (in prep.). In this region also the denso dwarfing 
gene was mapped [39], which may be identical with the semi- 
dwarf gene sdwl [40]. The second largest effect, after QHei.S42IL 
3Hb, was associated with QHei.S42IL-4H in S42IL-121. This 
QTL corresponds to QHei.S42IL-4H.a in Schmalenbach et al. 
[15]. In both studies the Hsp allele increased plant height by 18%. 
March et al. (in prep.) mapped a QTL for height for S42IL-121 as 
well. A third QTL (QHei.S42II^7H) with an increasing effect of 
the Hsp allele was detected on chromosome 7H. Here an effect 
that was already found in the studies of von Korff et al. [10], 
Schmalenbach et al. [15] and March et al. (in prep.) could be 
verified. Moreover two QTL where the Hsp allele had a decreasing 
effect on plant height [10] were verified. In S42IL-143 HEI was 
reduced by 11% (QHei.S42IL-lHb) and HEII (QHeii.S42IL- 
IHb) by 14%. Von Korff et al. [10] detected a QTL in the same 
region on chromosome IH. The flowering time gene HvFT3 is 
mapped in the same region and known to have a pleiotropic effect 
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on plant height [35]. QHei.S42IL-2H in S42IL-109 had reduced 
height by 9.5%. March et al. (in prep.) and Schmalenbach et al. 
[15] found the same effect in their studies. Von Korff et al. [10] 
found a similar effect in the region where the Hsp introgression of 
S42IL-109 was mapped. Moreover two candidate genes are 
mapped to the chromosomal region. These are the dwarfing gene 
sdw3 [41] and the flowering gene HvFT4, which is known to have 
an effect on plant height [35]. 

All of the HEI QTL in the present experiments were detected 
across treatments. Six out of eight QTL were also found in 
previous field and glasshouse studies. The QTL therefore seem to 
be very stable across locations as well as across treatments. 
Moreover they seem to be independent of the developmental 
stage. The present experiments, thus, allowed the verification of 
effects after six weeks that were previously screened in field 
experiments after flowering, indicating that phenot^ping juvenile 
plants may be predictive for adult plant performance, at least in 
regard to growth parameters. The high heritabUity of 76.8% 
supports this finding. 

For HEII two QTL coincided with previous studies. Besides 
QHeii.S42IL-lHb mentioned above, tiiis was QHeii.S42IL-4Hb 
where the Hsp allele reduced height by 14%. This QTL was also 
detected by von Korff et al. [10] and Wang et al. [35]. HeritabUity 
for digitally determined height was lower (61.4%) than for the 
manually measured one. Determining height by multiple mea- 
surements apparentiy was not an advantage here. However, this 
may change at a later stage of plant development. After shooting, 
the plant height is less subjected to bending of leaves and therefore 
can be measured more precisely by the imaging technique. 

Water Use Efficiency (WUE) 

Water use efficiency indicates how much biomass a plant can 
produce per unit water supplied. Thus, increased WUE has the 
potential to improve yield under drought stress conditions. 
Measuring WUE in regular greenhouse experiments is time- 
consuming. Therefore the high-throughput phenomics facility 
greatiy assisted in scoring of water use efficiency through 
automated watering of pots to specific weights. 

In this study the two S42ILs -117 and -129, with wild barley 
introgressions on chromosomes 4H and 6H, respectively, showed 
significant differences in WUE compared to Scarlett. Both ILs 
showed reduced water use efliciency compared to Scarlett. These 
ILs also produced less biomass. S42IL-117 and S42IL-129, thus, 
clearly carry unfavorable alleles for this trait. Chen et al. [42] 
pointed out that WUE itself is difficult to measure under field 
conditions and that a suitable tool to measure WUE efficiency is 
missing. Carbon isotope discrimination is a commonly used 
technique to measure WUE. Teulat et al. [43] used this method 
and identified QTL for WUE on chromosome 6H in a set of 167 
RILs from a cross between Tadmor and Er/Apm and likewise 
Diab et al. [44] identified a QTL for the same trait on 
chromosome 4H. The QTL detected in this study may correspond 
to the ones found in the studies mentioned before and suggest the 
results from both techniques are correlated. 

Compactness Integral, Shoot Area Top View Integral, Hull 
Area Integral (COM!, SATVI and HULI) 

The compactness of a plant describes how much of the hull area 
is covered by leaves. It was calculated as the ratio of SATVI to 
HULI. The more compact a plant is, the more ground cover it has 
with regard to the huU area. Two QTL were detected for this trait. 

SATVI and HULI showed a high correlation of r = 0.9, 
however, correlations between COMI and HULI and between 
COMI and SATVI were only moderate and negative. This 



indicates that in general, bigger plants take more space and have a 
lower compactness compared to smaller plants. In the present 
experiments this was observed by comparing drought stressed and 
well watered plants. Drought stressed plants showed on average a 
higher compactness than well watered plants. Jansen et al. [45] 
report the same effect on a study in Arabidopsis thaliana. 
Compactness shows negative correlation with all other traits 
evaluated in this experiment, with the exception of HEII (r = 0.16). 
An example for this is S42IL-117. This introgression line has a 
higher compactness, but reduced biomass, and other growth 
parameters compared to Scarlett. 

SATVI is one of the three parameters that control SAI and, 
thus, is highly correlated with this trait as well as with BMD and 
AGRI. As one may expect, for SATVI the same QTL were 
detected as for SAI. For HULI a total of six QTL were detected. 
Three of those may be due to high correlations in accordance with 
SATVI, AGRI and BMD. 

Caliper Length Integral (CALI) 

Caliper length describes the maximum diameter of the plant. 
For this trait six QTL were detected. Those were in accordance 
with QTL for HULI. This can be explained by the close 
connection of both traits. Hull area is taken as the basis to 
calculate caliper length and both traits are highly correlated 
(r = 0.94). CALI also shows positive and high correlations with 
AGRI, SAI, SATVI, HEI, and BMD. Plants with a large diameter 
cover a larger area, tend to be bigger, have a higher growth rate 
and a higher biomass than plants with a smaller diameter. 
Therefore, a lot of information on plant structure can be deduced 
from the plant diameter. 

Stress Indices 

Simple stress indices were calculated for each trait as the ratio of 
the mean plant performance under drought stress versus well 
watered treatments. In this study no QTL for a SSI was detected. 
Additionally the authors used two more complex stress indices 
(modified after Fischer and Maurer [46]), but were not able to 
detect QTL with those either. A stress index states how well a 
genotype performs under stress conditions relative to its perfor- 
mance under control conditions. Therefore, to see differences 
between genotypes for a stress index, a line by treatment 
interaction is necessary. If all genotypes show a similar growth 
reaction under stress and control conditions, the initially existing 
differences between the genotypes may be drastically reduced. In 
the present experiments line by treatment interactions were not 
significant and autocorrelations were high between the treatments. 
This may be the reason why no significant effect for the stress 
indices was found. This notion is supported by Wang et al. [47] . In 
their study on "mathematicaUy-derived traits in QTL mapping" 
the authors pointed out that an increased complexity of the genetic 
architecture of derived traits (e.g. stress indices) may reduce the 
power of QTL detection. 

High-throughput Phenotyping using The Plant 
Accelerator 

Determination of biomass by manual harvest is tedious and 
time-consuming. In addition, destructive harvest makes repeated 
measurements on the same plant impossible. Visual light imaging 
technologies applied in this study can solve these problems by 
utilizing the strong correlation between the projected shoot area 
and the actual biomass [27]. Imaging technologies ha\-c been 
successfully used in several studies in Arabidopsis thaliana, e.g. 
Granier et al. [48] and Leister et al. [49]. The first study 
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investigated nine accessions under different le\ els of water deficit 
in the phenotyping facility "PHENOPSIS". Reaction to water 
deficits was, amongst other traits, characterized by leaf area 
growth determined through images. The authors pointed out the 
importance of the automated watering in their experiment, which 
enables equal conditions for all plants. A characteristic that was 
also found very important in the present experiments. Leister et al. 
[49] described a first approach of using an image based technology 
for high-throughput growth analysis. They calculated plant area 
from top view images and found high correlations to plant fresh 
weight. 

In this study the sum of three two-dimensional pictures was used 
as a measure of plant biomass. In these experiments correlation 
between SAI and BMD and between AGRI and BMD were very 
high (r = 0.98). The results with six-week old barley plants proved 
that the sum of three pictures accounts sufficiently for overlapping 
leaves during early development. Rajendran et al. [30] found the 
same for T. monococcum. However, as Munns et al. [27] pointed out, 
accuracy may decline when plants become larger and produce 
multiple shoots. The results of the present study approved that 
The Plant Accelerator is suited to enable detailed growth analysis 
of barley plants. The prediction of biomass by the image-based leaf 
sum gave accurate results when comparing to actual biomass. 
Growth curves can give detailed information on differences in 
development of genotypes. For instance, the maximum of the 
absolute growth rate gives insight into the change from vegetative 
to generative phase of plant development. The present experi- 
ments ran only for six weeks. Therefore not all plants have reached 
this point. In future experiments this factor should be accounted 
for by adjusting the duration of the experiment. Automated 
imaging and the appropriate analysis pipehne make the detection 
of different developmental stages of plants feasible in high- 
throughput. With this technique it is possible to detect differences 
in stress responses between genotypes not only at different time 
points, but also to account for differences in development at those 
time points. 

Rajendran et al. [30] used non-destructive imaging to screen for 
different response mechanisms of T. monococcum to salt stress. In 
contrast to conventional Scdt stress experiments, where tolerance is 
measured as total biomass production of stressed plants compared 
to unstressed plants, the growth curves provided through daily 
imaging gave detailed insight into the tolerance mechanisms of the 
plants. While osmotically tolerant plants showed a constant growth 
rate, the growth rate of sodium excluders first dropped than 
increased after a couple of days. Moreover plant color was 
analyzed. No stress symptoms occurred on leaves in the present 
experiments. Due to littie variation between the genotypes no 
QTL was detected for leave color. However, in the experiments by 
Rajendran et al. [30], color analysis was successfully used to screen 
leaf damages due to high salt concentrations. 

Fluorescence imaging gives information on the health statues of 
a plant. It allows for detection of leaf senescence and necrosis. 
However, such symptoms were not observed in the present 
experiments and, thus, this parameter was not applied. Nonethe- 
less, the technique is readily available. In addition, near infrared 
(NIR) and infrared (IR) imaging may be useful for future plant 
growth evaluations and QTL studies. NIR enables the observation 
of the water status of a plant, while IR is used to determine shoot 
temperature. 



Conclusion 

In this study the use of a non-destructive high-throughput 
phenotyping platform was implemented to map QTL controlling 
\cgctative drought stress responses in barley. Several QTL where 
the exotic; H.sp allele had a positive effect on trait performance were 
detected. In particular, introgression line S42IL-121 showed 
improved growth under drought stress compared to the recurrent 
parent Scarlett. The line showed the same behavior in previous 
field experiments. Thus, this introgression line might be interesting 
for further breeding. 

Moreover, several QTL were detected where the Hsp allele had 
a decreasing effect on trait performance. Especially two QTL for 
water use efficiency might be interesting for further investigation. 
In future, interesting effects of S42IL-121 and other S42ILs will be 
fine mapped with already available high-resolution progeny [18] 
to further narrow down the QTL region and, ultimately, clone the 
underlying genes, which caused the observed QTL effects. 
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